Boundary quantum criticality 
in models of magnetic impurities coupled to bosonic baths 
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We investigate quantum impurity problems, where a local magnetic moment is coupled to the 
spin density of a bosonic environment, leading to bosonic versions of the standard Kondo and 
Anderson impurity models. In a physical situation, these bosonic environments can correspond 
either to deconfined spinons in certain classes of Z2 frustrated antiferromagnets, or to particles in a 
multicomponent Bose gase (in which case the spin degree of freedom is attributed to hyperfine levels). 
Using renormalization group techniques, we establish that our impurity models, which feature an 
exchange interaction analogous to Kondo impurities in Fermi liquids, allow the flow towards a stable 
strong-coupling state. Since the low-energy bosons live around a single point in momentum space, 
and there is no Fermi surface, an impurity quantum phase transition occurs at intermediate coupling, 
separating screened and unscreened phases. This behavior is qualitatively different from previously 
studied spin-isotropic variants of the spin-boson model, which display stable intermediate-coupling 
fixed points and no screening. 

PACS numbers: 75.20.Hr,74.70.-b 



I. INTRODUCTION 

The investigation of quantum phase transitions in 
strongly correlated electronic systems has been an active 
field of research during the last two decades. 1 Recently, 
there has been growing interest in so-called impurity 
quantum phase transitions - these correspond to qual- 
itative changes in the ground-state properties of discrete 
degrees of freedom (a spin 1/2 in the simplest case) upon 
changing their couplings to an external environment. On 
the one hand, impurity quantum phase transitions can 
be realized, e.g., in mesoscopic systems, such as quantum 
dot devices, which provide the high tunability required to 
access the transition point. On the other hand, the study 
of impurity problems is motivated through the dynamical 
mean- field theory 2,3 (DMFT), which maps lattice models 
of strongly correlated electrons - assuming a momentum- 
independent self energy - onto effective self-consistent 
impurity models. In this language, certain bulk quan- 
tum phase transition correspond to impurity transitions 
in Anderson or Kondo- type models. 

Impurity quantum phase transition occur in different 
contexts and setups. One possibility is a single mag- 
netic impurity coupled to an unconventional fermionic 
host. A paradigmatic example being the so-called pseu- 
dogap Kondo model, 4 originally motivated by the study 
of magnetic impurities in d-wave superconductors, where 
a spin 1/2 is coupled to the spin density of fermionic 
quasiparticles with a power-law density of states (DOS), 
p(uj) oc \ll>\ t ■ The lack of low-energy states (as com- 
pared to the standard Kondo model 5 ) causes a compe- 
tition between the local-moment dynamics and the for- 
mation of a strong-coupling ground state and leads to 
rich quantum critical behavior. 4,6-8 Other fermionic im- 
purity models display quantum phase transitions due to 
the competition between differently screened states, with 



the two-channel Kondo 9,10 and the two- impurity Kondo 
models 11-13 being popular examples. The majority of 
these models were recently reviewed in Ref. 14. 

Quantum impurities with a bosonic environment con- 
stitute a distinct class of models which can feature quan- 
tum phase transitions as well. A well-studied model is 
the spin-boson model, originally introduced to capture 
dissipation effects at the quantum level. Here, a generic 
two-level system (dubbed spin) is subject to the com- 
peting influence of a transverse field (causing tunneling 
between the two levels) and a longitudinally-coupled set 
of harmonic oscillators (providing friction). In the case 
of ohmic dissipation, this results in a Kosterlitz-Thouless 
phase transition at zero temperature, tuned by a varia- 
tion of the dissipation strength. 15 Note that the ohmic 
spin-boson model is equivalent to a Kondo model with 
anisotropic spin exchange terms, and hence does not con- 
stitute a distinct universality class. (This is different for 
sub-ohmic damping, see Refs. 14,16.) The topic of impu- 
rities in quantum magnets has motivated studies of the 
so-called Bose-Kondo model (which may also be viewed 
as a spin-isotropic spin-boson model with zero transverse 
field) - this model has a a stable intermediate-coupling 
fixed point and no screening; 17-21 this will be reviewed 
briefly in Sec. III. We note that the above-mentioned 
bosonic impurity models can be extended to include a 
fermionic environment as well, resulting in Bose-Fcrmi 
Kondo Hamiltonians, where the interactions of the impu- 
rity with the bosonic and fermionic baths compete. 22-26 
Models of this type have been proposed to describe quan- 
tum criticality in heavy-fermion compounds within ex- 
tended DMFT, 22 and could as well be realized in certain 
quantum-dot setups. 27 

We have recently introduced 28 a distinct class of 
SU(2)-symmetric bosonic impurity models, where the en- 
vironment consists of spin-1/2 objects, i.e., spinons con- 
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stituting the elementary excitations of certain frustrated 
quantum magnets. In contrast to the spin-boson type 
models mentioned above, where the environment always 
plays the role of a fluctuating (classical) magnetic field, 
our models display true quantum dynamics and allow 
singlet formation accompanied by complete screening. 

The purpose of this paper is to give a comprehensive 
renormalization group (RG) analysis of various versions 
of these bosonic Kondo models, thus providing details 
which were left out in the short report of Ref. 28: (i) 
The first example concerns the elementary bosonic ex- 
citations (spinons) in the Z2 spin-liquid state of a tri- 
angular Heisenberg model. 29 . This model will be shown 
shown to display an impurity quantum phase transition, 
bearing some resemblance to the fcrmionic pscudogap 
Kondo model, (ii) A second example that follows the 
same philosophy concerns the setup of a "magnetic im- 
purity" embedded in an optical lattice filled with a mul- 
ticomponent bosonic gas. 30,31 Here, the bosons are not 
collective modes originating from an underlying magnetic 
state, but are rather real bosonic atoms. In both cases, 
we derive the relevant quantum-field theory and study it 
using RG combined with epsilon expansion techniques. 
Although the structure of the theories is reminiscent of 
fermionic Kondo models, there are important differences: 
we demonstrate the importance of potential-scattering ef- 
fects at the impurity location as well as the influence of 
bulk interactions. 

Note that this paper will only study the RG flow at 
weak coupling; for situations with runaway flow to strong 
coupling - indicative of true screening - different tech- 
niques need to be applied. We have demonstrated this 
in Ref. 28, using strong-coupling toy models and large- 
TV techniques; this will not be repeated here. Further- 
more, we restrict our attention to situations where the 
bulk state does not break a local symmetry, i.e. we focus 
on paramagnetic phases including the quantum critical 
points. 

The plan of the paper is as follows. In Sec. II we 
summarize the order-parameter description of nearly- 
critical antiferromagnets, distinguishing collincar and 
non-collinear ordering. In the non-collinear case of in- 
terest to us, deconfined bosonic spinons 29 emerge which 
correspond to fractional spin 1/2 excitations in a Z2 spin 
liquid phase. 32-34 In Sec. Ill we review the results of 
earlier work 1 on a magnetic moment in a collincar 
antifcrromagnet, described by a conventional <f> theory 
with triplon excitations, leading to a spin-isotropic spin- 
boson-like model without screening. Sec. IV is devoted 
to the central topic of this paper, namely the physics of a 
magnetic impurity in a Z2 spin liquid. We present a mi- 
croscopic derivation of the impurity action, which we an- 
alyze by RG methods, extending our earlier calculation 28 
to two-loop level. Furthermore we propose a way to prop- 
erly incorporate the bulk interactions using a mapping to 
an effective bosonic Anderson impurity model. Finally, 
in Sec. V we study how an impurity two-level system 
couples to an atomic gas of canonical bosons carrying an 



internal degree of freedom (pseudospin). Again, we de- 
termine the phase diagram and critical properties using 
RG. A brief outlook will close the paper. Details of the 
RG calculations are deferred to the appendices. There, 
we also point out that the RG flow derived in Sec. V 
is identical to the one of a fully asymmetric fermionic 
Kondo model, which may arise in the DMFT context. 



II. FIELD THEORIES OF NEARLY-CRITICAL 
QUANTUM ANTIFERROMAGNETS 

Different types of magnetic Mott insulators naturally 
call for different types of order-parameter theories to 
describe the physics near the transition towards the 
symmetry-broken phase. In general, introducing an or- 
der parameter (f> for the magnetic fluctuations near the 
ordering wavevector K, a gradient expansion of the ac- 
tion naturally leads to a standard multicomponent (j) 4 
theory (or a non-linear sigma model). However, care is 
required for non-collinear spin correlations, as shown be- 
low. 

Assuming a magnetic order characterized by a single 
ordering wave vector K, the local spin operators can be 
parametrized as 

Si = <Re [$e iR - Si ^ = n x cos(K ■ x{] + n 2 sm(K ■ x, t ) , (1) 

where Xi denotes the position of lattice site i, and 
(j)(x, t) = Hi + in 2 is a complex vector field in spin space 
which describes the order-parameter direction and varies 
slowly in space and time. 



A. Collinear magnetic order 

For a state with static collincar order we have fixed 
vectors ni j2 with Hi x n 2 = 0; the (Si) in Eq. (1) on 
all sites i point parallel or antiparallel w.r.t. a common 
axis. The undoped insulator La2Cu04 is of this type (if 
we neglect the small spin canting due to Dzyaloshinski- 
Moriya interactions), with in-plane ordering wave- vector 
K = (tt,tt) (the lattice spacing is unity here and in the 
following) . For this specific value of K on a square lattice 
the order parameter is independent of n 2 . 

Near quantum criticality, one is lead to the usual 0(3) 
quantum 4 theory in imaginary time: 



S b = J d d xdr |a^| 2 + |9 r ^| 2 + 50 (^ 



I2\2 



(2) 



For K = (71", 7r) the vector field (f> is now real; for commen- 
surate K, cj) is complex, but the phase degree of freedom 
is gapped. This theory has been studied at length in 
the literature, 1 ' 29 and its critical properties are relevant 
for the quantum phase transition from a Neel state to a 
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gapped paramagnet in several classes of unfrustrated an- 
tifcrromagnets such as the bilayer square-lattice Heisen- 
berg model and coupled spin-ladder models. 

For general incommensurate K, one has to parametrize 
<f> = ne* e , with a real vector n and a phase field O. 
Interesting critical behavior may obtain here 35 , which is 
not in the focus of this paper. 

B. Non-collinear magnetic order 

Non-collinear order is present if n\ x n 2 ^ in (1); the 
order may be coplanar or fully three-dimensional. The 
former situation is e.g. realized in triangular quantum 
antiferromagnets, where the (Si) now lie in a plane in 
spin space, rather than along a single axis. 

The simplest case of a non-collinear ordered state is 
the one where 

ni • n 2 = ; n\ — n 2 — 1, (3) 

corresponding to a spin spiral. An elegant way to resolve 
the constraints is via a paramctrization of the order pa- 
rameter according to 32 

4> = ni +in 2 = e aal z a ,^-z T , (4) 

with the single constraint \z a \? = 1 (here a =f, J,), a 
is the vector of Pauli matrices, and e CTO -/ is the fully an- 
tisymmetric matrix. The objects z„ represent complex 
numbers (and hence are not canonical bosons), they de- 
scribe a mapping from the four-dimensional unit sphere 
S3 to S2. Obviously, this representation is unique up to 
a local Z2 transformation z a — ► r\ z a with 77 = ±1. 

Thus, in the non-collinear case described by (3,4), the 
long-distance fluctuations near a magnetic quantum crit- 
ical point are described by a quantum field theory in 
terms of bosonic spinons: 

S = - j d d xdr{\d x z a \ 2 + \d T z a \ 2 ) 

+ f d d xdTK(x,T)(Y,\ Z °\ 2 - 1 ) (5) 

where the integration over the auxiliary field re enforces 
the constraint on \ z <r\ 2 - I n Eq. (5) an additional Z2 
gauge field, corresponding to the above-mentioned gauge 
redundancy, has been omitted - this is justified provided 
that the gauge field is in a deconfined phase with massive 
fluctuations, which is the case for the triangular antiferro- 
magnet. Note that the paramagnetic phase of (5) is a Z 2 
spin liquid with gapped, deconfined bosonic spinons. 32-34 
Let us note here that the fluctuations in a collinear 
magnet may also be written in terms of bosonic spinons, 
leading to a Schwinger boson representation: 36 ' 37 

S{x) = Y J ^)^Ax), (6) 

a. a' 



with the constraint J^a \ z a\ 2 = 1- Note the difference 
with Eq. (4): in contrast to the non-collinear case, the 
accompanying gauge redundancy here is of U(l) symme- 
try. In such a situation, the gauge-field fluctuations are 
massless and cannot be ignored in a proper description 
of the bulk magnetism. The associated compact gauge 
theory is in general confining, 38 implying that bosonic 
spinons are bound in pairs and leading to the standard 
<j> 4 theory (2) for triplon fluctuations. 39 

U(l) spin liquids with deconfined spinons may never- 
theless arise: a special situation here are so-called de- 
confined critical points, 40 where a description in terms 
of bosonic spinons coupled to a non-compact U(l) 
gauge field applies (which implies the suppression of 
hedgehogs 41 ). Impurities characterized by a gauge charge 
and embedded in a U(l) spin liquid have been studied in 
Ref. 42; the interplay with a magnetic impurity moment 
would be interesting to study, but is beyond the scope of 
the present paper. 



III. A MAGNETIC IMPURITY IN A 
QUANTUM ANTIFERROMAGNET WITH 
COLLINEAR ORDER 

Magnetic impurities embedded in a nearly critical 
collinear quantum antifcrromagnet were analyzed in de- 
tail in Rcfs. 17-21. The bulk action Sb (2), describing 
the physics close to the transition between the ordered 
Neel state and the disordered spin-gap state, is supple- 
mented by the coupling to an impurity quantum spin 1/2 
5, leading to: 

Simp ~\~ ^Berry [S] , 

5 imp = 70 / drS{T) ■ ${x = 0, r) (7) 

where <fi(x, t) is the real order-parameter field of (2), and 
<SBcrry[S1 encodes the dynamics of the impurity spin S. 
A weak-coupling RG analysis 17 ' 18 leads to the RG beta 
function for the dimensionless impurity coupling 7: 

/3(7) = dwi = "y +7 + Til" (8) 

where e = 3 — d, fi is a rcnormalization scale, and the 
sign of the beta function is such that negative terms are 
RG relevant in the infrared. In Eq. (8), g denotes the 
renormalized bulk coupling in <S b (2) which flows as in the 
standard <j> 4 theory. 43 From Eq. (8) one finds a fixed-point 
value of the impurity coupling as 7* 2 = e/2 + C(e 2 ). The 
critical fluctuations of the bulk plus impurity problem 
are thus controlled by a fixed point with couplings g* 
(bulk) and 7* 2 (impurity), which arc both of order e, i.e., 
perturbatively accessible near d = 3 bulk dimensions. 

Most importantly, the flow of 7 (8) near 7* is infrared 
stable (!), in sharp contrast to (for example) the im- 
purity quantum critical point of the pseudogap Kondo 
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model. 4 ' 7 ' 8 The absence of runaway flow in Eq. (8) sug- 
gests the absence of a strong-coupling phase with com- 
plete screening, even at the bulk critical point with gap- 
less bosonic modes. In a sense, this is not unexpected, as 
Eq. (7) describes the classical fluctuations of an effective 
magnetic field coupled to the impurity, being unable to 
quench the entropy. 

As a consequence of (8), no quantum phase transi- 
tions occur in the model (7). The impurity properties are 
controlled by a stable intermediate-coupling fixed point, 
which describes a non-trivially fluctuating fractional-spin 
state. Near this fixed point, the impurity susceptibility 
and entropy obey 18 



HmTximp 



Ci(d), lim SUp = C 2 {d) . 



T^O 



(9) 



with universal constants C\(d), 62(d) (which depend on 
dimensionality d and impurity spin size S only). We 
note that the problem of an impurity in a collincar mag- 
net can equivalently be formulated using a non-linear 
sigma model for the bulk 19 - this leads to an expansion 
in e = d — 1, with similar impurity properties. Thus, the 
fractional-spin state is present for all 1 < d < 3 - for 
d = 2 the field-theoretic predictions 18,19 have been veri- 
fied by extensive Quantum Monte Carlo simulations. 20,21 
Paranthetically, we note that the fixed point at 7* is 
unstable w.r.t. breaking of the underlying SU(2) symme- 
try: a stable intermediate-coupling fixed point still exists 
in the XY case, whereas flow to strong coupling occurs 
in the Ising situation. 24,25 



IV. A MAGNETIC IMPURITY IN A Z 2 SPIN 
LIQUID 

As above, we intend to couple a magnetic impurity to 
the local order-parameter field of the bulk, i.e., by a term 
jo J drS (t)-(P(x — 0,t), but now the low-energy dynamics 
of <j)(x, t) will be represented in terms of spinon fields z 
introduced in Sec. II B. 



A. Low-energy theory 

It is easy to see that an appropriate choice of a linear 
combinations of the spinons defined in Eq. (4), namely 



4 = ^( z T+* z l)> 



(10) 



enables us to write expression (4) in an alternative, but 
very convenient way: 

4>{x = Q)=^t{n 1 +in 2 ) = z':^-zy (11) 
Note that this rotation leaves the bulk action invariant. 



The appropriate model to describe a magnetic impu- 
rity in a topological spin liquid with a Z2 gauge structure 
is thus 

S = ~ g j d d xdT(\d xZr7 \ 2 + \d T z a \ 2 ) 

+ j ' d d xdTK(x,T)(^\ Za \ 2 -l) 



(12) 



+ jo J drS-z;^-z a ,(x = 0) 

+ <5Berry[S1 ) 



and was proposed recently by us. 28 

The most important feature of model (12) is that 
the impurity is coupled to a bilinear of the elementary 
bulk excitations (i.e., the spinons). This has to be con- 
trasted to the impurity theory with collinear bulk or- 
dering, where the impurity spin couples linearly to the 
magnetic (triplon) modes. The consequences of this dif- 
ference will be examined below and turn out to be cru- 
cial for the fate of the impurity near the bulk critical 
point. Thus, apart from the strong constraint on the bulk 
spinons, the model (12) strongly resembles the fermionic 
Kondo model and constitutes its straightforward exten- 
sion to a bosonic environment. 

In the following, we concentrate on the impurity be- 
havior at the bulk critical point, i.e., for gapless spinons. 
In the gapped case, weak impurity coupling will be irrel- 
evant, and the impurity will decouple from the bulk in 
the low-energy limit. However, in analogy to fermionic 
Kondo models with a gapped density of states and broken 
particle-hole symmetry, 14 a first-order quantum transi- 
tion will take place upon increasing jo, beyond which a 
strong-coupling state is realized. 



B. Weak-coupling RG analysis 

A RG treatment of the action presented in Eq. (12) 
can be performed perturbatively in j - In addition, the 
constraint z* a z a — 1 is relaxed and treated on the 
mean-field level (this is equivalent to giving the spinons a 
mass shift). The influence of such bulk interaction terms 
will be considered in Sec. IV C. At zero temperature, the 
renormalized spinon mass vanishes at criticality, and the 
spinons become soft. This is the regime we focus on, and 
which is described by the simpler action: 

S = ~ g J d d xdr(\d xZa \ 2 + \d T z«\ 2 ) 



+ jo drS-z*^-zAx = Q) 



+ S: 



Berry 



[S]. 



(13) 



In order to use diagrammatic methods it is convenient to 
represent the impurity spin in terms of Schwinger bosons 
(one could equally well choose to work with Abrikosov 
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fcrmions). Restricting ourselves to an impurity spin 1/2, 
we introduce thus bosonic operators b a which obey the 
following constraint 



(14) 



This allows for a faithful impurity spin representation 
according to 



(15) 



In terms of these bosonic operators, we can formulate the 
full impurity action: 44 

S = 1 J d d xdr z* a (-flg - d 2 T ) z a + J drb a (d T + A) b a 

+ j J dTb a ^b a ,z* a ^-z (x = 0) 

+ v j. J drz^x^O), (16) 

where we implement the constraint (14) exactly by taking 
the limit A -> oo 45 - 46 

In Eq. (16), we have included a potential scattering 
term vq on the impurity location x = 0. In contrast to the 
fcrmionic case, such a term is generated perturbatively in 
the RG flow. This implies that for finite j, v = cannot 
be a fixed point of the RG - this is intimately related to 
the bosonic nature of the bulk, as bosons cannot show 
particle-hole symmetry. 

The RG treatment starts by determining the tree-level 
scaling dimensions of the couplings, yielding: 



dim[jo] = dim[w ] = d — 2 = e . 



(17) 



Both couplings are marginal in two dimensions; this 
allows for a controlled RG procedure, which we 
shall perform here using the standard field-theoretic 
scheme. Renormalized couplings are introduced as jo = 



s§^' u e ^T-^ an< ^ v ° = g§d' Lt e ^lr t '' wnere 3d denotes 
Sd = r ^| ( 2 4i)J/r At two-loop order (see App. A) we 

find the following beta functions 



0U) = z3+vj + Y + 0(f), 



(18) 



where e = d — 2 for non-interacting bulk bosons. (This 
corrects a sign error in the result published in Ref. 28.) 
It is important to note that f3(j) [/3(f)] does not contain 
even (odd) powers of j, respectively. This is the con- 
sequence of the invariance of the action (16) w.r.t. the 
transformation 



Jo 



-jo and z a 



(19) 



which also implies that ferromagnetic and antiferromag- 
netic impurity coupling arc equivalent. (Note again here 



that the z are not canonical bosons - the corresponding 
operators simply commute.) 

The solution of the two-loop flow equations is de- 
picted in Fig. 1. The weak-coupling RG flow has four 
fixed points: the trivial (Gaussian) fixed point LM with 
{ v *ii*) — (0,0); a potential scattering fixed point PS, 
given by (v*,j*) = (— 2e,0); and a pair of critical fixed 
point, located at (v*,j*) = (— e — ^,±-^=) and denoted 
QCP. 

Let us quickly discuss the physical implications of the 
present RG flow (but note that certain changes will arise 
from the consideration of bulk interactions, see below) . 

For e > and j ^ 0, the flow indicates two sta- 
ble phases, namely the local-moment phase (LAI), and 
a strong-coupling phase at large \j\ (the Hamiltonian is 
symmetric under jo — > —jo) and large negative v. The 
two phases are separated by a continuous quantum phase 
transition, controlled by QCP, with a correlation length 
exponent given by 



1 

v 



(20) 



characterizing the vanishing crossover energy scale near 
criticality. The pure potential scattering problem, j = 0, 
has two phases as well, which can be easily understood 
from the exact local spinon Green function: 



G(iv) 
where Go{iv) 



1 + (V4)GoM 
d d k g 



I 



(27T) d v 2 + k 2 



(21) 
(22) 



Values v > — 2e flow to zero (the denominator in Eq. (21) 
is regular), whereas v < —2e induces a bound state (a 
pole appears in the Green function of the bosons). The 
physical properties of LM are that of a decoupled impu- 
rity, i.e., susceptibility and entropy obey Ximp = 1/(4T), 
S'imp = In 2 in the low-temperature limit. The strong- 
coupling regime cannot be analyzed using the present 
RG; a combination of large- N methods and strong- 
coupling models indicates that true screening can occur 
(for both signs of jo), with a finite impurity susceptibil- 
ity and vanishing entropy. 28 Finally, at the critical point 
(QCP) properties similar to Eq. (9) obtain. 

For e < 0, the weak-coupling regime is unstable, and 
the flow generically is towards strong coupling - this can 
be related to the bare boson density of states, Eq. (22), 
diverging as v — > 0. Nevertheless, there is still a QCP, 
which separates the strong-coupling (screened) regime 
from a pure potential-scattering phase at positive v (con- 
trolled by a now stable PS fixed point). 

For € = no quantum phase transition occurs, and 
there is logarithmically slow flow to strong coupling (see 
below for corrections to this picture coming from bulk 
interactions). 

Overall, the RG flow is strikingly different from an 
impurity which is coupled to gapless spin-1 bosons 
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(a) 6 > 



(b)e<0 




Figure 1: RG flow of the spinon Kondo model (16) at the 
bulk critical point, (a) e > 0; (b) e < 0. (e = d — 2 for 
non-interacting spinons; note that e may be non-zero in d = 2 
owing to bulk spinon interactions.) QCP denotes the bound- 
ary quantum critical point and the blue line is the separa- 
trix, separating a strongly coupled phase from either a local 
moment (LM) regime or a potential scattering (PS) a zero 
Kondo coupling. The curves have been calculated from (18) 
for e = ±0.1. Note that the flow is symmetric w.r.t. j «-» —j. 



(Sec. Ill), where the coupling 7 flows towards a stable 
intermediate-coupling fixed point. 



C. The role of bulk interactions 

Within the above weak-coupling RG the bulk interac- 
tions (represented by the constraint \z a \ 2 = 1) were 
only treated at the mean- field level. Unfortunately, a 
consistent RG fully including the bulk interactions can- 
not be performed: the appropriate bulk non-linear sigma 
model is perturbatively accessible near d = 1 (while soft- 
ening the bulk interaction allows an expansion in 3 — d, 
see below) , but expansions around both d = 1 and d = 3 
are incompatible with the bare scaling dimension of the 
impurity coupling jo, which is marginal in d = 2. (This 
is in contrast to the situation in Sec. Ill where both bulk 
and boundary couplings were marginal at d = 3.) 

The treatment of bulk interactions in the present 
spinon case is, however, physically important, both at 
weak coupling (to be discussed below) and at strong cou- 
pling - this was discussed in some detail in Ref. 28, but 
requires further work which is beyond the scope of this 
paper. 

Let us qualitatively discuss the role of bulk spinon 
interactions at weak coupling. The scaling dimensions 
of j and v are determined by the local bulk spin cor- 
relations. In the interacting bulk theory, this decay is 
given as x( T ) x \/T d ~ 1+ri where 77 is the anomalous di- 
mension. (Note that r\ = 1 for non- interacting spinons, 
whereas r\ = in the case of a standard 4 theory con- 
trolled by a Gaussian fixed point.) This implies that 
dim[jo] = dim[v ] = {d~3 + i])/2, i.e., the e = d— 2 in the 
RG equations (18) should be replaced by e = (ef— 3+77) /2. 
The QCP in the incommensurate magnet has r\ > 1 in 



d = 2 (for details see Refs. 32,33), i.e., a boundary quan- 
tum phase transition as shown in Fig. 1(a) is indeed pos- 
sible for the physical two-dimensional case. 

We note that this analysis (i.e. including r\ in the 
scaling dimensions of j and v) is in principle correct 
at one-loop order; at higher orders interference pro- 
cesses between bulk and boundary interactions can be 
expected. No systematic treatment of those is possi- 
ble in the present formulation. Therefore, we propose 
a more rigorous way to incorporate bulk interactions, 
which leads us to a bosonic version of the Anderson im- 
purity model. 



1. Derivation of the bosonic Anderson model 

Our goal is to find a representation of the model (16) 
which allows to perturbatively control both bulk and 
boundary interactions at the same time. To this end, 
we reformulate the impurity part, i.e., switch from a 
Kondo to an Anderson impurity model - such a scheme 
has proven extremely successful for the pseudogap Kondo 
model. 

Two transformations will convert the original model 
(16) into a new one, where bulk and boundary interac- 
tions have a common critical dimension: (a) an inverse 
Schrieffer- Wolff transformation (in analogy to the pseu- 
dogap Kondo case) , which transforms the Kondo into an 
Anderson model and renders the hybridization between 
impurity and bulk marginal in d — 3; (b) replacing the 
hard constraint on the spinons by a "soft" self-interaction 
of the type mo(|z,t| 2 — l) 2 _ this term is marginal in d = 3 
as well. (As will become clear, step (a) requires the exis- 
tence of a strong-coupling singlet state with screening in 
the original model.) With steps (a) and (b) we are led to 
propose the following field theory: 



- (\d x z a \ 2 + \d T z a \ 2 ) + u Q z*z tT z*,z (7 , 
9 



S = J d d xdr 

+ / dr \b (d T + A + e ) b + b a (d T + A) b a 



dr \b a boz a + h.c] 



(23) 



We have left out a bulk mass term oc z 2 , as we are in- 
terested in the bulk critical point. The impurity is now 
a three-level system (similar to an Anderson model with 
infinite repulsion), represented by three auxiliary parti- 
cles b a (a =T)i) and bo - note that bo may be inter- 
preted as the empty impurity state. The impurity tran- 
sition is tuned by the value of to'- f° r eo > the im- 
purity is in a doublet state (i.e. unscreened) whereas 
for eo < the state is a singlet (these statements ap- 
ply for wo — > 0, finite wo will shift the impurity critical 
point). The Hilbert space dimension is enforced by the 

constraint Q = J^a ^t^a + ^0^0 = Technically, this is 
again implemented via a chemical potential A — > 00 



45.4(i 
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The action (23) incorporates three main ingredients: 
(i) it contains an interaction term u$ in the bulk, which 
is marginal in d = 3; (ii) the hybridization wq is marginal 
in d = 3 as well; (iii) in the "Kondo" limit, i.e., for a 
large positive singlet energy eo, a Schrieffer- Wolff trans- 
formation allows to integrate out the bo field, and leads 
to the previous bosonic Kondo interaction (16). The 
model (23) at d = 3 is thus strongly reminiscent of the 
pseudogap Anderson impurity model (with infinite impu- 
rity Coulomb energy) and linear density of states. 8 (An- 
other bosonic Anderson model was recently introduced 
by Lee and Bulla, 47 which, however, displays different 
physics, as the model of Ref. 47 is spinless, but with fi- 
nite on-site repulsion.) 

We will proceed with deriving the RG flow for the An- 
derson model (23). 



2. RG analysis of the bosonic Anderson model 

The RG flow for both the bulk and impurity inter- 
actions, uq and wq, can be perturbatively controlled in 
e = 3 -rp . This allows a consistent calculation of the crit- 
ical properties of the fully interacting bosonic impurity 
model. The details of this calculation are presented in 
App. B. 

We introduce renormalized couplings according to 

wo = u% /-| — ; Zm w and uo — M 2£ -o^-?w u - Plac- 

ing ourselves at both the bulk and boundary critical 
points (!), we find the following flow equations to one- 
loop order: 

3 

(3(w) — ew — -^w 3 , 

(3{u) = 2eu- 12m 2 . (24) 

These RG equations describe the non-trivial flow of both 
bulk and boundary interactions at the critical point in 
d < 3. Higher-order contributions (not shown here) con- 
tain the above-mentioned interference processes and in- 
duce a feedback of the bulk interaction u into the impu- 
rity flow. 

We refrain from a more detailed analysis at this point: 
qualitatively, the RG approaches in Sees. IV B and IV C 2 
yield similar results; a quantitative comparison is diffi- 
cult as the two expansions are around a lower critical 
and upper critical dimension of the impurity problem, 
respectively. Let us, however, emphasize that we have 
demonstrated how a systematic expansion, analogous to 
the one in Ref. 18 for impurity models in collinear mag- 
nets (i.e. treating bulk and boundary interactions on 
equal footing) can be set up in the non-collinear case as 
well. 



V. A MAGNETIC IMPURITY IN A 
TWO-COMPONENT CANONICAL BOSE GAS 

In this last part, we briefly introduce and analyze the 
problem of a two-state system 30 (described by a magnetic 
pseudospin impurity) coupled to a two-component non- 
relativistic Bose gas. 31 This extends the kind of bosonic 
impurity models, introduced above in the context of 
quantum antifcrromagncts, to the field of cold atomic 
gases. 



A. The model 

Technically, the main difference to the situation in the 
previous section is that the bulk now consists of canonical 
bosons. The action assumes the form: 

S = J d d xdTa a (x, t) (<9 t — <9 2 + to) a a 
+ jo J drS ■ a a ^-a al (x = 0) 
+ ^ J dr a a a a (x = 0) 
+ S Bmy [S] . ( 25 ) 

where we consider canonical bosonic particles, a a , with 
two internal degrees of freedom (a =|,| here) and a 
mass to (i.e., chemical potential). In the field of ultra- 
cold trapped gases, usually performed with fixed particle 
number, to — > in the low-temperature limit. Interac- 
tions among bulk bosons are neglected; to lowest order, 
they lead to a temperature-dependent renormalization of 
to. The impurity couples to the bosonic spin density at 
site x = with a Kondo-type exchange interaction, but 
an Anderson-like Hamiltonian similar to Eq. (23) could 
be considered as well. 



B. Weak-coupling analysis 

In contrast to spinon-based models, we realize that 
there is no invariance under the transformation jo - * 
—jo, i-e., ferromagnetic and antiferromagnetic impurity 
coupling arc no longer equivalent. However, as in the 
spinon case, a potential scattering term is generated dur- 
ing the flow even for v = 0. 

In the remainder, we focus on the weak-coupling anal- 
ysis in the massless case, to = 0. (In the massive case, 
to > at T = 0, we again expect a first-order transi- 
tion as function of jo-) A tree- level dimensional analysis 
yields 

dim[j ] = dim[u ] = d — 2 = e, (26) 

similar to the model in Sec. IV B. We thus introduce 
dimensionless couplings according to jo = fJ'~ e 4~% L j an d 
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vq = /i e ~^%^ v and nn d to one- loop order 

m = £ i-y + f > 

/?(«) = e^ + ^ + ^j 2 . (27) 

Remarkably, this result is exact to all orders in perturba- 
tion theory (this was previously shown for non-relativistic 
4 -theory with local interactions 1,48 ). Technically, the 
reason for that lies in the diagrammatic structure of the 
problem, where only ladder-like diagrams can be differ- 
ent from zero in the limit T = 0. Higher-order terms 
are cancelled by the appropriate counterterms. Physi- 
cally, the exactness of Eq. (27) is tied to the fact that the 
expansion is performed around a trivial vacuum state. 

The RG flow, derived from Eq. (27), is depicted in 
Fig. 2. For positive (antiferromagnetic) j, it is qualita- 
tively similar to the one of the spinon model in Sec. IV B. 
From Eq. (27) one finds four fixed points: the local- 
moment (LM) fixed point at (v*,j*) — (0,0), a potential- 
scattering (PS) fixed point at (v*,j*) — (— 4e, 0), and two 
critical fixed points at (v*,j*) = (— e,e) and (v*,j*) = 
(— 3e, — e). For both QCP, the correlation length expo- 
nent is 

- = N , (28) 
v 

which again is exact to all orders in perturbation theory. 

The structure of the phase diagram for both posi- 
tive and negative e is thus similar to the spinon case in 
Sec. IV B. However, ferromagnetic and antiferromagnetic 
jo are no longer equivalent, suggesting the existence of 
two distinct strong-coupling phases for j < and j > 
(where true screening is likely to occur only in the latter 
case). 

For e = and j > 0, logarithmic flow to strong 
coupling occurs, with a screening temperature given by 
InT* cx — 1/j'o as in the conventional Kondo problem. 

Thus, the situation of a two-level impurity coupled to a 
two-component Bose gas, with a pseudospin-pseudospin 
coupling, allows for three distinct phases (given jo =/= 0) : 
a weakly coupled (unscreened) phase and two strong- 
coupling phases. The present analysis can be easily gen- 
eralized to an arbitrary number of pseudospin compo- 
nents or larger effective spin sizes; at weak coupling no 
qualitative changes are expected. However, the precise 
nature of the ground state at strong coupling depends on 
microscopic parameters of the Bose gas, such as spin size 
and the value of the repulsion between bosonic atoms. 28 
More elaborate strong-coupling methods are required to 
address this issue. 

Last not least, let us point out another remarkable 
feature of the above model (25): its RG flow is equivalent 
to a fermionic Kondo problem with full band asymmetry, 
i.e, with a chemical potential tuned to the band edge. 
This was considered previously, and we comment on 
this point in App. C. 




Figure 2: RG flow of the Kondo model with massless canonical 
bosons, Eq. (25). (a) e > 0; (b) e < 0. (e = d — 2 for non- 
interacting bosons; again, e may be non-zero in d = 2 due 
to bulk interactions.) The curves have been calculated from 
(27) for e = ±0.1. 



VI. CONCLUSION 

In this paper, we have introduced a new class of quan- 
tum impurity models, describing the dynamics of a two- 
level system (spin) coupled to an external bosonic spin- 
carrying environment. Two different situations have been 
studied in the weak-coupling regime, namely a spin 1/2 
magnetic impurity in a class of Zi spin liquids, and a 
local defect in a two-component Bose gas. 

For massless bulk particles, we have derived the renor- 
malization group flow of the impurity coupling constants. 
While some differences in the RG flow appear between 
the above two realizations, the following generic features 
appear: (I) an impurity quantum phase transition is pos- 
sible which separates a local-moment regime (decoupled 
impurity) from a strongly coupled impurity; (II) in con- 
trast to a magnetic impurity coupled to a collinear bulk 
magnet described by a standard 4 theory, there is no 
stable intermediate-coupling fixed point, but instead run- 
away flow to strong coupling. (For gapped bulk parti- 
cles, a first-order transition will generically appear, in 
analogy to gapped fermionic Kondo models.) At strong 
coupling, the nature of the ground state depends on mi- 
croscopic details (and was not investigated further here), 
but quenching of the impurity entropy may obtain. 28 

In summary, our results underline that the dynamics 
of impurity degrees of freedom depends crucially on the 
nature of the elementary excitations of the surround- 
ing environment. Hence, impurity properties may be 



9 



used as probes for exotic bulk excitations. On the ex- 
perimental side, CS2CUCI4 has been proposed to be a 
realization of a spin-1/2 frustrated antiferromagnet on 
an anisotropic triangular lattice, which possibly features 
spinon excitations. 50 ' 51 Studying the physics of dilute 
magnetic impurities in this material, e.g. the temper- 
ature dependence of NMR Knight shifts or of the impu- 
rity contribution to the susceptibility, will clearly help to 
resolve the nature of the bulk magnetic state. 



a. Weak-coupling RG procedure 



The tree-level scaling dimensions of the general cou- 
pling constants are 



dim[a ] = dim[/? ] = dim[7o] = d - 2 = e . (A2) 
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Appendix A: WEAK-COUPLING RG FOR THE 
BOSONIC KONDO MODEL 



Based on this analysis we introduce dimensionless cou- 
plings according to 



a 

A) 

7o 



2 Z a 
god Z b 

-e 2 z 7 
gb d z h 



(A3) 



In this appendix we derive the RG equations for the 
bosonic Kondo model with bulk spinons (the derivation 
is kept in a form which is also appropriate to treat the 
case of canonical bosons as well as electrons; the only 
occurring modifications manifest themselves in the val- 
ues of the corresponding Feynman diagrams). The RG 
equations of bosons are very different form those of elec- 
trons, which is the reason why we present a rather de- 
tailed analysis. The crucial difference resides in the fact 
that during the renormalization process a potential scat- 
tering term is generated. In order to properly deal with 
this complication we start with a generalized version of 
the model proposed in Eq. (12) which allows for fully 
spin-anisotropic couplings. This enables us to deal with 
the Kondo interaction and the potential scattering term 
at a time without worrying about newly generated terms. 

The action we analyze has the following very general 
form 

S = \ [d d xdTj24(-d 2 T -dt)z„ 

+ zZ f drl - & + A ) b ° 

+ a Q / dr ^2 z Z z v'bcr'ba 
+ A) / dT^z* a z a b a b a 

J cr 

+ 7o / dr ^2 z*zj> a ,b^ ■ ( A1 ) 



where /x denotes an arbitrary but fixed energy scale (the 
renormalization scale) and Sd = 1 where 0<j = FTdJTj 
with T(x) being the Eulcr function. We furthermore state 
that during the course of renormalization the Schwingcr 
boson field parametrizing the spin is renormalized ac- 
cording to 



b R = Z h - X b . (A4) 



We do not have to introduce renormalized fields for the 
bulk degrees of freedom since the influence of a single 
impurity on the bulk is an effect of the order 0(1/N), 
which does not renormalize bulk degrees of freedom. The 
prefactor can trivially be taken care of by rescaling 
the fields in Eq. (12) according to 




which has the only advantage of getting rid of bother- 
some factors in the course of the calculation. The gen- 
eralized renormalized action (with appropriate countert- 
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fi a 



l)a 



1)0 



1) 7 



Figure 3: Graphical representation of all the occurring ver- 
tices; a is denoted as a black circle, its counterterm is repre- 
sented by a big black circle in a box; the small grey dot and 
the big grey dot with a cross represent /3 and its counterterm; 
the small (big) grey box (grey box with a cross) represents 7 
(the counterterm of 7) . 



erms) reads: 

Sr = t / ddxdT £ z ° (~ d * ~ z ° 

J a 

+ /i~ e 7 J dr ^2 Zo,aZo,ab a <b a < 

+ J dTj2~b a (d T + \)b a + J2j dr{Z b - l%d T b a 
+ fi~ e a(Z a - 1) J dr ^2 Zo,a Z 0,<T'ba'ba 
+ n~ € f3(Zf3 - 1) / dT^ZQ^zo^bgba 

J a 

" e l( z i ~ 1) J dr ^2 z *),<j z Q,<?b al b a , . 

(A6) 



This is the appropriate starting point to iteratively con- 
struct the counterterms by demanding the action to be 
finite at a certain arbitrary but fixed renormalization 
point. The graphical representation of the vertices and 
the counterterms is shown in Fig. 3. The vertex a is de- 
noted as a small black dot whereas the associated coun- 
terterm [Z a — 1) a is represented as a big black dot in 
a box. The interaction vertex (3 and its counterterm 
(Zp — 1) (3 are represented as a small grey dot and a big 
grey dot with a cross in the interior, whereas 7 and the 
associated counterterm (Z 1 — 1)7 are graphically repre- 
sented as small grey box and a big grey box with a cross 
inside, respectively. In the diagrammatic expansion this 
enforces us to consider 6 interaction vertices for a consis- 
tent treatment. 

The next step is to specify the renormalization condi- 
tions which will give the theory a finite limit once the 



cutoff is sent to infinity according to 
T\{iv-\ = iv = <S) = 



aw 



= 1 



IV — fl 



r 4 , , 

88 88 



fi a , 
At" 6 /?, 
^~ e 7 > 



(A7) 



where the subscript denotes the appropriate renormaliza- 
tion point, which we do not specify here. The first line 
fixes the impurity propagator, whereas the second line 
fixes the vertex functions. We once again stress the fact 
that the bulk properties will not receive singular correc- 
tions in any order of perturbation theory. 



b. Result to one-loop order 

To lowest loop order we can concentrate on the renor- 
malization of the vertex functions since there is no prop- 
agator renormalization at one loop level (they will play 
a role at two loop level). To lowest loop order there are 
no propagator renormalizations. We choose a graphical 
representation of the renormalization conditions fixed in 
Eq. (A7) . The renormalization conditions for the vertex 
functions are shown in Fig. 4, Fig. 5, and Fig. 6. Struc- 
turally there are only two types of diagrams occurring, 
whose most divergent part in an expansion in 1/e yields 



dx- 



x(l + x) 



-fi 



1 



O(e 



(A8) 



By dressing the integral values with the appropriate in- 
teraction vertices we can determine the Z factors to low- 
est order as 



Z\ 



3 



(2/3 + 2 7 ) - , 

a 2 \ 1 
— + 27 - . 

7 ) e 



(A9) 



c. Result to two-loop order 

In order to derive the Z factors to two-loop order we 
only calculate the numerical values of the topologically 
distinct diagrams shown in Fig. 7. The actual vertex cor- 
rections are obtained by properly dressing the numerical 
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Figure 4: Graphical representation of the renormalization 
condition shown in Eq. (A7) determining the one-loop coun- 
terterm (Z a — Xy- 1 ' (big black dot in box). 




Figure 5: Graphical representation of the renormalization 
condition shown in Eq. (A7) determining the one-loop coun- 
terterm (Zp — l)' 1 ^ (big grey dot with a cross). 



value of the integrals with the corresponding interaction 
vertices. All possible diagrams contributing to the two- 
loop renormalization factors are shown in Figs. 7 and 8. 
The diagrams contributing to Z 7 are not shown explicitly, 
since they can be obtained from the diagrams contribut- 




Figure 6: Graphical representation of the renormalization 
condition shown in Eq. (A7) determining the one-loop coun- 
terterm [Z 1 — l)' 1 ' (big grey box with a cross). 



ing to Zp by replacing grey dots by grey boxes. We start 
with the trivial diagrams (a) and (6) which are just the 
squares of the one-loop diagrams 

(a) = (b)=^(±+O(e )\ . (A10) 

The diagrams (c)...(j) are given by 

(c~j)=M 2e I™ dx ^*^ [°°dy-^, (All) 
Jo 1 + x J 1 + y 

which to leading order yields 

(c.... J )=M 2e (^+<^°)). (A12) 

The remaining diagram (k) evaluates to 

(*) = M 2e j°° dxdy- — T~T77 T \ 
Jo [x + y+l)(x + y) 

= -M 2e (^ + C(e )). (A13) 

We will re-encounter the same integral within the calcula- 
tion of the field renormalization factor. In order to arrive 
at the final two loop Z-factors we dress the integral val- 
ues with the appropriate interaction vertices. After some 
calculation we find 




+ |(2(7 2 +« 2 ) + ^) ■ (A14) 

The next step consists of the propagator renormalization, 
whose recursion relation is shown in Fig. 9. We calculate 
the diagram shown in Fig. 9 (a). The dependence on the 
external frequency can be extracted by calculating 

E(w„ - A) - S(0) 

= {iv n - A)/i 2e / dxdy- — ■ — — — ■ — - 

Jo (x + y+l){x + y) 

= -{ iVn -\)^(L +O {e )j . (A15) 

If we dress this integral with the appropriate vertices 
and apply the renormalization condition shown in Fig. 9 
(note that iV n = iv n — A) we find 

Z b = l + ^{a 2 +l3 2 + 1 2 ) . (A16) 
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Figure 7: Graphical recursion relation for the two-loop contribution to the renormalization factor (Z a — l)' 2 ' (last term). Note 
that the counterterms appearing as internal vertices in the diagrams take their one- loop values (Z a - l) w a, (Zg - 
[Z 1 — 1) 7- ( a )-(k) denote topologically distinct diagrams. 




Figure 8: Graphical recursion relation for the two-loop contribution to the renormalization factor (Zp — l)' 2 ' (last term). Note 
that the counterterms appearing as internal vertices in the diagrams take their one- loop values (Z a - l) (1) a, (Z - l) (1) /3, 

(Z-y ~ 1) (1) 7- 



b) c) In order to derive the RG beta functions we use the fact 

that the bare couplings have no knowledge about the 
" - x *" " > renormalization condition, which can be expressed as 



da dp d<y 
H—=H—=H— = 0. (A17) 
d/i an djU 



Figure 9: Recursion relation for the impurity field renormal- 

ization factor (Z b - 1) (crossed dot). This leadg tQ & get q£ coupled equations 
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— ea + a- 



da 



91n(#0 91n(#* 
^/3(a) + a h^W) + « H^^M = 



5/3 



#7 



9 In 



e/3 + /3 



£7 + 7 



9a 
9a 



/3(a) +0 



dln(^) 91n(f^ 

,J -m+P } // /3( 7 )=0, 



9/3 



97 



/3(a) +7- 



9 In 



6) 



9/3 



/?(/?) + 7" 



9 In 



(t) 



97 



/3( 7 ) = 



(A18) 



where the introduced /3 functions correspond to /3(a) = 
H-j-a etc. These coupled equations can be solved itcra- 
tively, leading to 

/3(a) = ea + a (a 2 + (/? - 7 ) 2 ) + 2a(/3 + 7 ) , 



/3(/3) - e/ ? + 2/3^— +/?J + /?(^a 2 - 
0(7) = e 7 +(a 2 + 2 7 2 )-a 2 (/?- 7 ) 



a z 7 



(A19) 



In order to arrive at the actual bosonic Kondo model 
(Eq. (12)) with potential scattering we have to make the 
following identifications 

a=|, 0=^ + 1, and 7 = + J . (A20) 

The flow equations of the physical couplings jj_ , j z , v are 
related to the intermediate quantities a, /3, 7 according 
to 



- 20(a) , 
/3( Jz ) = 2[/3(/3)-/3( 7 )] 
/?(«) = 2[/3(/3)+/3( 7 )] 



(A21) 



with the interaction parameters chosen as shown in 
Eq. (A20). This leads to the following flow equations 



ej_L + vj± + 



3±3z , J± 



<(<■) = + £ + £ + 



(A22) 



or in the spin-isotropic situation 



ej + + y 



(A23) 



which are the RG equations (18) quoted in the main text. 

The RG equations for the model of canonical bosons, 
Eq. (27) in Sec. V, are obtained from the same diagram- 
matic expansion (albeit with different numerical values 



for the diagrams). The observation that the RG equa- 
tions are exact to all orders comes can be easily under- 
stood by noting that the counterterms cancel all higher- 
order contributions because the only non-vanishing dia- 
grams are ladder-like (in electronic language one would 
call these diagrams the Cooper diagrams), see Ref. 1. 



Appendix B: RG FOR THE BOSONIC 
ANDERSON MODEL 

Within this section we derive the RG equations for 
the bosonic version of the infinite-?/ Anderson model 
of Sec. IV CI. After a Fourier transformation for the 
quadratic bulk part, the action of the model (23) reads 

+ u J d d xdTz'*z,jz'*,z al 

J dTb (d r + \)b + J dTb a (d T + X)b a 

, J dr[b a b a z a (T)+h.c] . (Bl) 



+ 
+ w 



The flow equations for the bulk interaction w and the hy- 
bridization wo will not mix to one-loop order. In the fol- 
lowing we will ignore the renormalization of the mass dif- 
ference between the different impurity levels [eo in (23)], 
i.e., we work at criticality. We introduce renormalized 
couplings according to 



wo = 



uo = A* 



2c 



4 Z u 

u , 



and 



S d g 2 Zl 



(B2) 



where e = ^=^. It is interesting to note that technically 
we have a version of the RG of a model which is at its 
upper critical dimension (both for the bulk and the im- 
purity), whereas the RG for the bosonic Kondo model 
is reminiscent of an RG at the lower critical dimension. 
The vertices of the theory are shown in Fig. 10. 
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(a) 



(b) 



(c) 



Figure 10: (a) Interaction vertex between bulk and impurity 
in the bosonic Anderson model; the wiggly line denotes the 
particle 60, the dashed line is b a and the full line stands for z a 
(incoming arrow) or z * (outgoing arrow), (b) and (c) denote 
self-energy corrections to the local levels b a (b) and 60 (c). 



The calculation for the renormalization of the hy- 
bridization w completely parallels the analysis in the 
pseudogap infinite-C/ Anderson model, which was pre- 
sented in Ref. 7, and therefore without showing further 
calculational details leads to the same flow equation, 
namely 



dlV 3 n 

— = ew — —w 

dm/i 2 



(B3) 



In addition, we have to calculate the perturbative correc- 
tions to the bulk interaction vertex, which is identical to 
the usual 0(4)-</> 4 theory. In our convention the result is: 



du 
din \l 



-2eu+12u 2 



(B4) 



The hidden 0(4) symmetry of the theory becomes ap- 
parent upon writing the bosonic spinons as a 4-vector 



Appendix C: FERMIONIC KONDO PROBLEM AT 
THE BAND EDGE 

Remarkably, the two-level system coupled to the two- 
component bosonic gas studied in Sec. V turns out to 



follow a weak-coupling RG flow indentical to that of a 
fermionic Kondo problem, where the chemical potential 
is tuned to the (e.g. lower) band edge (i.e. with full 
particle- hole asymmetry). The equivalence can best be 
understood in terms of the local density of states - the 
only quantity which determines the impurity behavior - 
which in both cases (canonical bosons and fermions) is 
of the form p(ui) = p \uj\ r 8(uj). 

The strict particle-hole asymmetry of the electrons 
implies that there are only particle excitations. Intro- 
ducing dimensionless couplings according to standard 
practice 7 we can derive the following flow equations for 
the fermionic problem: 



m 



= rj 



3 



2 



3 n 1 2 

4 4 



(CI) 



These equations are contained in the earlier analysis of 
Ref. 49, and are identical to those obtained for the prob- 
lem of canonical bosons in Sec. V. Again, the equations 
are exact to all orders in perturbation theory (!). 

For a general fermionic metallic (r = 0) Kondo prob- 
lem with band asymmetry (i.e. the chemical potential 
not being pinned to the center of the band) one has to 
perform a two-step RG procedure: Initially, one has to 
follow the fully asymmetric scaling, Eq.(Cl). Once the 
remainder of the band is symmetric, the usual poor man's 
scaling equations take over - here the potential scattering 
term is a marginal operator. 

We note that the fully asymmetric fermionic Kondo 
problem can be relevant in the context of DMFT for 
Hubbard models: The transition from the undoped to 
the doped Mott insulator upon variation of the chemical 
potential happens precisely when the chemical potential 
reaches the gap edge (provided the transition is continu- 
ous). 
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